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The major thrust of the computational analysis of turbomachinery to date 
has been the steady-state solution of Isolated blades using mass-averaged Inlet 
and exit conditions. Unsteady flows differ from the steady solution due to 
Interaction of pressure waves and wakes between blade rows. To predict the 
actual complex flow conditions one must look at the time accurate solution of 
the entire turbomachine. 


ANALYSIS 

Many of the numerical tools used In the analysis of Isolated blades can 
be used for time-accurate analysis. Here the quasl-three-dlmenslonal Euler and 
thin- layer Navler-Stokes equations are solved for unsteady turbomachinery 
flows. These equations are written In general coordinates for an axlsymmetrlc 
stream surface, and they account for the effects of blade- row rotation, radius 
change, and stream- surface thickness (ref. 1). 

A four-stage Runge-Kutta scheme based on the work of Jameson (ref. 2) Is 
used to predict time-accurate results. Body-fitted C-type grids were used In 
this work and were generated using the GRAPE code (GRIds about Airfoils using 
Poisson's Equation) developed by Sorenson (ref. 3). 

To predict the Interaction between a stator and a rotor In turbomachinery, 
one must pass flow Information to an Interface that acts as a moving boundary 
between the two computational grids. A nonconservative Interface formulation 
Is used In this procedure. The solution must be Integrated using a constant 
minimum time step based on the computational domain of two blades. The data 
management necessary to update the stator and rotor flowflelds will be dis- 
cussed. Currently the code Is limited to solving stator/rotor configurations 
that have equal pitch blades. 


RESULT S 

The solution procedure has been applied to two test problems. The first 
Is a cascade and the second Is a turbine rotor from the space shuttle main 
engine (SSME). Euler and Navler-Stokes results will be presented. 

A model turbine stage made up of two Identical cascades of NACA 0012 air- 
foils was used to develop and test the Interface and data management routines. 
Figure 1 shows the C-grlds used In the computation where the second blade row 
Is moving downward. A converged solution can be seen In the periodic loading 
diagram of the rotor (fig. 2). 
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The stator Inlet Mach number Is 1.84. Figure 3 shows relative Mach number 
contours when the blade rows are aligned. Detached bow shocks form ahead of 
the stators. Interact at mid-passage, and reflect obliquely back to the stator 
surface. A strong curved shock forms at the stator trailing edge, and the 
shock curvature generates an entropy layer that convects downstream into the 
rotor. The stator flowfield Is symmetric but the rotor flowfleld is asym 
metric due to the incidence of the relative flow. 

Figure 4 shows relative Mach number contours after the rotor has moved 1/4 
pitch downward. The flow between the stator and rotor is subsonic so that 
pressure waves from the rotor can affect the stator trai 1 ing- edge shock, which 
has moved upstream. The entropy layer from this shock produces a large asym 
metry In the rotor flowfield. lhe flow reaccelerates to supersonic in the 
rotor passage and produces strong oblique shocks off the rotor trailing edge. 

Previous steady-state calculations of the SSMF first-stage fuel- side tur- 
bine blade (ref. 1) showed a reverse flow region on the pressure surface 
(fig. 5). lime- accurate results are being generated to determine how that 
separation is affected by the upstream stator wakes, lhe grids used for this 
calculation are shown in figure 6 where the second blade row is the rotor of 
Interest. Euler results are shown in figure 7. There Is very little interac- 
tion between the blade rows. Viscous results will be shown in the presenta 
tion. Currently the numerical procedure works on blade row configurations of 
equal pitch. The actual SSME stator/rotor blade count Is 41:63, which is close 
to the 2:3 count shown In figure 8. A future goal of this work is to use this 
solution procedure to solve problems of unequal blade count like the one shown 
here. 
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Figure 7. - Static pressure 




